
restore,'./ancillary/grs_kda_data.sav',/ver

omni_lbd2rz,grs.glon,grs.glat,grs.dis*1.d3,r,z

plothist,z,bin=10,x,y,thick=3

yf = mpfitpeak(x,y,A,nterms=3)

xp = findgen(4001)/10. - 200.

cgOplot,xp,gauss_1(xp,A),thick=3,color='red'

al_legend,/top,/left,box=0,['FWHM = '+string(A[2]*2.355,format="(F0.1)")+' pc']

END
